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Abstract 



We study the modifications of synchrotron radiation of charges in a stor- 
age ring as they are cooled. The pair correlation lengths between the charges 
are manifest in the synchrotron radiation and coherence effects exist for wave- 
lengths longer than the coherence lengths between the charges. Therefore the 
synchrotron radiation can be used as a diagnostic tool to determine the state 
(gas, liquid, crystal) of the charged plasma in the storage ring. We show 
also that the total power of the synchrotron radiation is enormously reduced 
for crystallized beams. This opens the possibility of accelerating particles to 
ultra-relativistic energies using small-sized cyclic accelerators. 
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I. INTRODUCTION 



Ion-beam crystallization is an exciting and relatively new field of physics in which a new 
state of matter is sought for. Namely, ions which rapidly circulate in a storage ring and 
are cooled are expected to form geometrically-ordered structures (crystals) which have a 
density much smaller than normal crystalline solids []I]|§. Although great effort is currently 
invested in achieving such crystals , and cooling techniques were significantly improved 
||, there is still no clear-cut experimental evidence for them. It is hoped, however, that 
crystalline beams will be produced in the near future. 

Synchrotron radiation, on the other hand, is a very well-established field of physics 
that has been investigated continuously from the early days of particle accelerators. Many 
synchrotron sources are operating around the world (e.g., DESY (Hamburg, Germany), 
NSLS (Brookhaven, USA), KEK (Tsukuba, Japan)), and many applications already exist 

!• 

It is the purpose of this paper to establish a link between beam crystallization and 
synchrotron radiation. This link is two-fold: 

• To use synchrotron radiation and modifications thereof in order to detect the creation 
and existence of beam crystals. This is required since for fast beams direct detection 
methods are difficult to implement 0. Thus synchrotron radiation can be used as an 
indirect diagnostic method to detect the formation of beam crystals. The diagnostic 
methods discussed below are also applicable to liquid and gaseous beams. 

• Even more importantly, once beam crystals are formed, they can be used to modify 
the synchrotron radiation with respect to the ordinary incoherent case. Therefore one 
can achieve dramatic suppression and enhancement effects of synchrotron radiation 
using crystallized beams 0. In particular, the total power that is radiated from an 
equi-spaced circulating chain of particles is much smaller (in the appropriate limits 
[[]]) than the radiation from the same number of randomly-located particles. This 
opens the possibility for accelerating particles to ultra-relativistic energies with little 
radiation loss, which is currently the main limitation of circular electron accelerators. 
Thus the suppression of synchrotron radiation by beam crystallization may eventually 
lead to the construction of smaller-sized circular electron accelerators. 

In the following we shall detail the connection between beam crystals and synchrotron 
radiation. It is important to emphasize that currently researchers are trying to obtain 
beam crystals of heavy ions that can be cooled with electrons and lasers. For heavy ions, 
however, the synchrotron radiation is small. Even for protons, e.g., the lightest of the "heavy 
ions", the synchrotron radiation is about a factor 10~ 13 smaller than for electrons with the 
same energy. Thus, we expect that realistically the effects predicted in this paper will be 
important only for liquid or crystallized electron beams. This, however, poses the challenge 
of obtaining crystallized electron beams. Thus we hope that the ideas put forward in this 
paper will motivate experimentalists to work towards obtaining crystallized electron beams. 
In any case, however, the analysis presented below applies to any species of charged particles. 
Thus the theory can in principle be verified for ion-beam crystals. Also we stress from the 
outset that the effects discussed below go beyond what is known as "coherent synchrotron 
radiation" which is the coherent enhancement of synchrotron radiation of small electron 
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bunches for wavelengths that are longer than the bunch size (see, e.g., [pHilf and Sec. [IV 
below) . 

The paper is structured in the following way. In Sec. [TT] we discuss the general theory 
that underlies the suppression and enhancement effects of synchrotron radiation. In Sec. 
I IT] the theory is applied to the three phases of a coasting charged-particle beam that occur 
in practice: Gaseous, liquid and crystal. In Sec. |IV| the necessary modifications of the 
theory for a bunched beam are discussed. In Sec. [V] we present analytical and numerical 
results concerning the suppression of synchrotron radiation by a crystalline beam. Finite 
temperature effects are discussed explicitly. In Sec. |VT| we discuss our results and conclude 
the paper with proposals for experimental applications of the effects discussed in this paper. 



II. GENERAL THEORY 



We consider N charged particles with charge q circulating in a circular storage ring 
of radius p with velocity v. The charges are assumed to be coherent with a reference 
circulating charge, but are allowed to have constant time lags Atj as well as constant spatial 
displacements Afj from the reference orbit. According to the theory of radiation of moving 
sources, the total power that is emitted from the N charges is given by [|l2|,|n| 



(2.1) 



n=l 



where is the power that is emitted with frequency cu n = nto = nv/p due to a single 
circulating particle, and g n is the form factor of the beam. The explicit expression for 1^ 
is 
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where (3 : 
functions 
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v/c, c is the speed of light, 7 
The form factor is given by 



1/vT — (3 2 , and J n are the ordinary Bessel 



N 



3=1 



N 



exp[m(0 i - <f) f ) 



(2.3) 



and the angles 4>j are given by JTJj 



uAtj H — n ■ Ar 
P 



(2.4) 



where n is the unit vector pointing from the center of the ring to the observation point. We 
observe that in the form factor the role of the time delays and the spatial displacements 
is equivalent. Thus we restrict ourselves hereafter to time delays only. This simplifies the 
calculations and gives qualitatively the same results. It is also compatible with the current 
experimental trend according to which linear ion crystals (one-dimensional crystals) are 
sought for. We shall denote in the following the phase differences by 9j = uAtj. 
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Suppose now that we treat the quantities 9{ as random variables distributed according 
to the normalized probability density P{9\, . . . , 9^). Then the expectation value of the total 
power is 

oo 

(I {N) } = J2(9n)li 1 \ (2.5) 

n=l 

where 

,2tt N 



(g n ) = J d9 1 ■ ■ ■ d6 N P(9 1 , ...,9 N ) £ exp[in(^ - 9 r )}. 
i,.i'=i 



(2.6) 



3,3 

This can be rewritten as 

•-2tt 



(g n )= [ ^ dAe inA i? 2 (A), (2.7) 

J-2tt 



where 

R 2 (A) = J d9i---d9 N P(9 h ...,9 N ) £ 5[A -(%-%,)] (2-8) 

is the two-point correlation function, i.e. the (non-normalized) chance of finding a pair of 
#'s a distance A apart. Therefore we conclude that the crucial quantity that determines 
the modifications of synchrotron radiation due to coherence effects is R2, our main object 
of study. The physics of the particle beam (temperature, structure) is reflected in R2 and is 
consequently linked to modifications of the synchrotron radiation. 

Before applying the above formulas, we make some further simple manipulations. To 
avoid complications with the 27r-periodicity we define 

_R 2 (A) = 2[R 2 (A) + R 2 (2n - A)]. (2.9) 

Using the relation R 2 (A) = R 2 (—A), easily derived from ( [2.8|) , we obtain 

tg\ = TdA cos(nA)£ 2 (A). (2.10) 
Jo 

This can finally be recast as 

(g \ = N+ r dA cos(nA) J R 2 (A), (2.11) 
Jo 

where _R 2 (A), defined in [0, n], is the two-point correlator that does not include the "diago- 
nal" part 2NS(A), emerging from the j = j' terms of R 2 . 

An important special case is the case of independent particles, i.e. 

N 

P(9 1 ,...,9 N ) = Y[P 1 (9 j ), (2.12) 

where Pi is the (normalized) one-point density of the particles. For this case the resulting 
form factor is 



(g n ) = N + N(N-l) 



2 77 



d9P 1 {9)e 
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(2.13) 
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III. APPLICATION TO COOLED PARTICLE BEAMS 



In the following we shall study a few representative situations of a particle beam as 
it is being cooled and crystallized. We shall qualitatively infer the form of the two-point 
correlator R 2 for each of the cases, and calculate the resulting form factor of the synchrotron 
radiation. Thus the focus in this section is on the spectral modifications of the synchrotron 
radiation expressed by the behavior of (g n ). We shall show that the modifications of (g n ) as 
the temperature is lowered defines an excellent tool for the diagnostics of the thermodynamic 
state of the beam. The modifications due to bunching are considered in Sec. |V[ The 
suppression of the total emitted power is discussed in Sec. |V|. 

We start with a very hot particle beam. In such a case we expect the particles to be 
completely independent. Therefore equations (2.12) and ( |2.13| ) apply. For a particle beam 
that fills the whole ring (coasting beam) we expect on the basis of symmetry a uniform 
distribution 



(3.1; 



The form factor becomes 



l ) = N. 



n 



1,2, 



(3.2) 



Thus, for a hot coasting beam, 



(JW) = NI™ 



(3.3) 



This is what we expect from totally incoherent radiation of N particles. 

For beams that are bunched, a typical shape is a Gaussian. The resulting density is 



V2 



TXO- 
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(3.4) 



where 8q is the location of the center of the bunch and a is its angular width. The summation 
over m is to ensure the 27r-periodicity. For 8q ^> a, 2tt — 8q ^> a only the m = component 
in ( |3.4j ) is significant. The resulting form factor is 



(g n )=N + N(N-l)exp(-n 2 a 2 ) 



(3.5) 



In the bunched case, therefore, in addition to the incoherent term N, we also obtain a term 
that represents the coherent synchrotron radiation for low harmonics n ^ l/er. One obtains 



qualitatively the same results for other shapes of the bunch |L7 , |I8 



The above results are well-known and form the basis of the field of "coherent synchrotron 
radiation" in which enhancement of the radiation is predicted [|TJ,|17],p,|rjJ and experimentally 
measured [p| [TTj] due to the collection of the charges (electrons) into small bunches. In our 
case this applies to the limiting case of a bunched but very hot beam of particles in which 
the particles within the bunch are uncorrelated. 

In the following we shall introduce the correlations between the particles as the beam is 
cooled and use the full expressions (|2.5|) - (|2.11 ) rather than fl2.13|) . These correlations are 
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neglected in the field of coherent synchrotron radiation since the beam is assumed to be very 
hot. But the correlations become more and more important as the beam is being cooled. To 
simplify the treatment, we focus in this section on coasting beams. In Sec. [IV] we introduce 
the necessary modifications to describe bunched beams. 

If the temperature of the beam is moderately high, we expect that the particles start to 
show repulsion from each other. That is, they will avoid the vicinity of each other due to the 
Coulomb repulsion and act like a non-ideal "gas" . This can be described phenomenologically 
by 



i?f S (A) = Cl R° 2 (A) 



1 — exp 



A 2 



2a 2 



(3.6) 



where 



Rr, 



(A) 



< A < ii 



7T 



(3.7) 



is the trivial two-point correlator for a uniformly coasting beam of independent particles 
and a is the (angular) "hard-core" scale of repulsion. Interpreting (|3.6| ), we modified R® by 
a narrow "dip" of width a near A = (the "correlation hole") such that i?2(0) = (total 
repulsion at A = 0). The constant c\ is for normalization. It is approximately (1 — a/ v^r) -1 
for small values of a (a <C tt). Actually, since we are still in the high temperature regime, we 
need to assume that a -C 2n/N = do, i.e. the hard core repulsion occurs on scales smaller 
than the mean distance between the particles. For N ^> 1, which is the interesting case 
here, the two conditions on a are consistent. We note that a depends on the temperature 
and increases as the temperature decreases. For a< 1 we obtain the following form factor: 



<sr> 



Ar N(N-l)a 
N v , ' exp 



'2tx 



n 



1,2,... . 



(3.8) 



Since the form factor is a non-negative quantity, we immediately infer an upper limit on a: 



< 



N 



2,71 



(3.9) 



This result is intuitively clear since the hard core cannot be larger than the mean distance 
of the particles. It is also compatible with the assumptions above concerning a. Physically, 
we observe that there is a suppression of the synchrotron radiation for n ^ 1/a (lower 
harmonics). We can therefore estimate the hard-core scale (and hence the temperature) 
from the coherent modifications of the synchrotron radiation for the gas-like state of the 
particle beam. We note that the overall suppression effect is small. For small n values, 
where the reduction of the emitted power is largest, the relative suppression with respect to 
the incoherent case amounts to only 



- N 



N 



27va 



<C 1. 



(3.10) 



The situation is illustrated in Fig. [I]. 
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FIG. 1. The two-point correlator (a) and the form factor (b) for the "gaseous" state of the 
beam, described by (3.6) and Q3,8|), respectively. We used the parameters a = dg/10,N = 100. 



When the temperature becomes smaller such that the Coulomb energy is comparable to 
the thermal energy, we expect the particle beam to become somewhat ordered and to form 
a liquid-like plasma. The partial order is a precursor to crystallization. In particular, the 
(angular) distance on which the repulsion between particles is manifest is dg, and the order 
effects should persist over a few mean distances. The two-point correlator is qualitatively 
given by 
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Rp(A) = c 2 R° 2 (A) 

where c 2 ~ (1 — 1/N)^ 1 is a normalization factor. As before, we assume N 1. The 
above two-point correlator displays a strong repulsion for small distances (/^(O) = 0) 
as well as oscillations that persist for a few mean distances. It eventually reaches the 
asymptotic limit of uncorrelated particles. Thus it represents an intermediate situation 
("liquid") between the slight mutual repulsion ("gas") treated above and long-range order 
("crystal") discussed below. The above expression was obtained as a result of an exact 
calculation for a one-dimensional chain of particles with logarithmic repulsion by Dyson in 
the context of Random Matrix Theory ||19|| . The resulting form factor is 

(g l ^)=mm(n,N). (3.12) 

Thus, the suppression effect is very prominent in this situation, and there is effectively 
complete suppression of the synchrotron radiation for small values of n (see also Fig. 0). In 
terms of wavelength, the suppression is felt for wavelengths that are comparable or longer 
than the mean distance between the particles. Comparing (|3.8[ ) and fl3.12| ) we conclude that 
as the order becomes more manifest (temperature decreases) the suppression effect becomes 
more prominent, but the onset of suppression is shifted to longer wavelengths. 
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FIG. 2. The two-point correlator (a) and the form factor (b) for the "liquid" state of the beam, 
described by (pID and ( gig ). 



As crystallization takes place, long-range order effects become important. We consider 
in the following the simplest crystal, namely the linear chain. To describe the situation 
we assume a distribution function that corresponds to a thermal distribution of small dis- 
placements around the crystalline state with only nearest-neighbor interactions taken into 
account for simplicity: 
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c 3 exp 



N 



-v ^2(^+1 -<Pj) 2 



3=1 



(3.13) 



Here C3 is a normalization constant, rj = (q 2 N 2 )/(167r 3 €odk-BT), T is the temperature and 
d = 2-np/N is the mean distance between the charges. The small displacements ipj are 
defined as follows: 



(j-l)d e , j = 2,3,...,N, 



The exponent rj can be rewritten as 



N 
2^ 



2 / typical potential energy 
I typical kinetc energy 



N 
2^ 



r. 



(3.14) 
(3.15) 
(3.16) 



(3.17) 



where T is the plasma parameter in one dimension f2| 



Aire^dk^T 



(3.18) 



We note that the assumption of only nearest-neighbor interactions is not severe since for 
small displacements the interaction with the n'th neighbor reduces as 1/n 3 . 

In this case it is easier to obtain the form factor directly, without explicitly calculating 
the two-point correlator. A straightforward but lengthy calculation gives the following result 
for the form factor 



(g^) = N + 2j:(N 
1=1 



1) cos 



27m/ 
N 



exp 



n 2 n 2 l(N - I) 
fiV3 



(3.19) 



In order to obtain the above result we assumed T ^> 1, i.e. a cold beam. This is a necessary 
condition for crystallization. In order to interpret this result, we consider two limiting cases. 
If the maximal exponent in (|3.19| ) (as a function of I), given by n 2 ir 2 / (4IW), is much smaller 
1, we can replace the exponential in ( |3.19|) with 1 and get 



than 



N 2 
, 



A^ divides n, 
otherwise. 



(3.20) 



This means that for very cold crystals there is a total suppression of the radiation for all 
harmonics, except the ones that are divisible by the number of particles N. For these special 
harmonics we get total constructive interference. The suppression of the leading harmonics 
results in an enormous reduction of the total power emitted by the synchrotron radiation 
(see Sec. 0). In case crystallized electron beams can be produced, this effect gives rise to the 
possibility of significantly reducing the synchrotron radiation, currently the main limitation 
for circular electron accelerators. We mention in passing that the result (|3.20|) can also be 
obtained directly from calculating the form factor for a completely frozen crystal [J7 , [12 , 13 . 
The other limit of ( |3.19| ) is for the first exponential factor (I = 1) to be already small, such 



that only the first term needs to be considered. That is, for n tt /(FN ) ^> 1, we obtain 
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<^ y > « N + 2iVcos (^) exp ("^) , (3-21) 

which describes small "ripples" over the incoherent radiation, with decaying amplitude that 
has oscillations with period N. In Fig. ||| we plot the numerically computed form factor 
( P-19| ) for the specific cases N = 10 6 , Y = 10 1 , 10 2 , 10 6 . For T = 10 6 we see a series of sharp 
peaks located at n/N — 1,2, ... . This is expected since in this case T is very large and thus 
( p.2U| ) holds approximately for the range of n shown in Fig. |3|. For smaller values of V we 
observe a transition from sharp peaks to decaying ripples. Even in the case V = 10 6 the 
sharp peaks will eventually die away. 
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FIG. 3. The form factor for the crystalline state of the beam, described by ( p,19| ). We considered 
the cases N = 10 6 ,T = 10\l0 2 ,10 6 . 

The depression of g c * y at n m 1 can be computed analytically. Expanding the exponential 
factor in ( |3.19| ) to first order in 1 /Y we obtain 



cry _ , 1 
— Oi ~ . 



(3.22) 



This is in perfect agreement with the results displayed in Fig. |3|. 

The above results concerning the crystalline state indicate that the plasma parameter Y 
can be determined from the form factor of the synchrotron radiation (provided iV is known). 
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This defines a useful diagnostic tool for measuring the temperature of the crystal. We also 
conclude that crystalline beams can be applied to selectively suppress and enhance harmonics 
of the radiation, achieving up to total suppression (g n = 0) or total constructive interference 

To summarize this section, we have shown that the synchrotron radiation and its modi- 
fications with respect to the incoherent state are strongly connected with the physical state 
of the beam. The form factor reflects the important scales and can be used to diagnose the 
state of the beam ( "gas" , "liquid" , "solid" ) as well as its temperature. 



IV. MODIFICATIONS FOR BUNCHED BEAMS 



Experimentally it is sometimes useful to work with bunched beams in which the particles 
occupy only a small fraction of the ring. Thus we consider in this section the modifications 
of the above theory for bunched beams. These modifications are straightforward. It turns 
out that only the lowest harmonics (up to n ~ 2tc/ (bunch angular length)) are affected. 
Qualitatively this can be understood by examining equation ( |2.6| ), since the bunching will 
be felt only for values of n such that n(9j — 6jr) < 2tt. This yields the above estimate. In 
the following we detail the theory quantitatively. 

We start with the gaseous phase and consider a narrow bunch of iV particles with an 
(effective) angular width a = 2n/Q,Q 3> 1. In order to be specific we shall assume that 
the shape of the bunch is a Gaussian, and that the (one-point) charge density is given by 
equation fl3.4| ) above. In the absence of correlations, the two-point correlation function of 
the Gaussian bunch reads 



^ B (A) 



N(N 



sf- 



Tia- 



+00 



(A + 2ixmf 
Ao~ 2 



(4.1) 



from which we calculate the form factor (|3.5| ). In order to include the hard-core repulsion 
between the charges, we operate as in the coasting case and modify R^ 3 with a narrow dip 



Dgas. bunch 
-^2 



(A) = c 4 i? 2 GB (A) 



1 — exp 



A 2 



2a 2 



(4.2) 



As before, C4 ~ 1 to leading order in N. When calculating the form factor, the first term in 
the brackets gives (|3.5|). For the second term, we can use R§ B ~ 1 since we assumed a<a. 
Hence, we obtain 



(9. 



gas, bunch 



) = N 



N 2 a 



exp 



+ N(N - 1) exp(-nV) 



(#f s > + N(N - 1) exp(-nV^ 



(4.3) 



That is, the form factor of the gaseous coasting beam contains an additional enhancement 
feature for low harmonics, n < Q. This is suggestive, because of the scale separation 
between the length of the bunch and the hard-core scale, a <C a. For the liquid phase a 
similar analysis applies. We need to replace the term R\ in ( p.ll[ ) with i? GB , and similar 
considerations will lead to the conclusion that we get the same type of enhancement of the 
low harmonics due to bunching 
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liq, bunch \ 



) = (^ q ) + iV(iV-l)exp(-nV 2 ) 



(4.4) 



For the crystalline state (linear chain) with finite temperature we model the bunch by 
adding two stationary charges at both ends of the bunch. These charges do not radiate and 
serve only for confinement. To make the calculations tractable, we assume only nearest- 
neighbor interactions. Lengthy but straightforward calculation yields the form factor 



N-l 

^cry.bunch^ = N + 2 ^ (N - I) COs(nldg) eX P 
1=1 



n 2 djl(N +1-1) 
4T(iV + 1) 



(4.5) 



For small values of n we replace the exponents by 1 and obtain 



(9n 



cry, bunch 



sin 



(nn/Q) 



[™/{QN)\ 



sin 



(4.6) 



This is the form factor of a frozen linear crystalline bunch. In particular, it exhibits an 
enhancement for the low harmonics n < Q. If n is so large that only the first term is 
significant, we essentially recover the result ( |3.21| ). Results for the case 
lOOv^ and T = 10 1 , 10 2 are shown in Fig. |j. The parameters were chosen such that dg is 
the same as for the coasting case. We observe that significant enhancement indeed occurs 
for the lower harmonics, which is essentially independent of the temperature as suggested 
by (|4.6|) . Otherwise, the form factor (normalized by the number of charges) is the same as 
for the coasting case. 
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FIG. 4. The form factor for the crystalline state of a bunched beam, described by Q4,5| ). We 
considered the cases N = 5V2 x 10 3 , Q = W0V2, T = 10 1 (upper plot), T = 10 2 (lower plot). 

To summarize this section, we investigated the modifications that result from the bunch- 
ing of the particle beam. In all cases we found that a significant enhancement occurs for 
the low harmonics n < Q. Otherwise we get qualitatively the same results as for a coasting 
beam. 



V. TOTAL POWER 



In Sees. [TTI| and [TV] we concentrated on a discussion of the form factor g n of the beam. 
We showed that important information on the thermodynamic state of the beam is already 
contained in g n . The total emitted power, however, the subject of this section, depends on 
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the interplay between g n and the partial power levels 1^ of a single radiating charge (see 
Eq. fl2.1p ). The total power 1^ of a single radiating charge is given by 

J« = £/« = -^/3V. (5.1) 

This result agrees with Larmor's well-known formula for the total radiated power of a single 
charge in the nonrelativistic limit We introduce the parameter s = f3j. It characterizes 
the three relativistic regimes important for the discussion in this paper: Nonrelativistic 
(s <C 1), relativistic (s ~ 1) and ultra-relativistic (s ^> 1). With the help of the total power 
Q5.1D we define the normalized power levels 

Since the purpose of this section is to discuss suppression effects in the total emitted 
synchrotron-radiation power, we define the suppression factor 

n=l 

In the case of iV incoherently radiating charges we have a(N,/3) — 1. A suppression ef- 
fect corresponds to a (AT, /3) < 1. Enhancement of synchrotron radiation corresponds to 
a(N,0) > 1. 

The behavior of 7^ as a function of n is the key for understanding the suppression effect 
of the total emitted synchrotron power. It is qualitatively different in the three relativistic 
regimes (see Fig. [|). For s C 1 we have (3 <C 1 and 1^' decays exponentially in n. This 
is illustrated in Fig. ||(a). It shows I$p as a function of n for s = 0.1. Expanding (|2.2|) to 
leading order in {3 we obtain 

3(n + l)n 2n+1 

(2n + l)(2n)! P ' P K ' 

We verify that ij 1 ^ ~ 1 in this limit. Using Stirling's formula we obtain 
/m~ 3(n+l)y^ fep\ 2n 

ln 2(2n + l)0F/^ ^ 2 J ' P<<i ' n>>i ' (5 - 5j 

which proves the exponential decay of 1^ for large n. The result (|5.5|) is also shown in 
Fig. |5](a). The exponential decay for large n persists in the case s ~ 1, albeit with a much 
smaller decay constant. This is illustrated in Fig. ^(b). In this case we also have an analytical 
approximation. It is given by ]14 



/ i \ 2 rt 

^"T^h , 1S7,»>7-. (5-6) 

The analytical approximation (|5.6|) is shown as the dashed line in Fig. |5](b). It describes the 
numerical data very well. The same figure also shows that a qualitative change with respect 
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to the nonrelativistic case (Fig. |5](a)) occurs only for small n where 1^' starts with a near- 
zero slope. In the ultra-relativistic case (s > 1) the behavior of 1$ changes qualitatively. 
For small n it shows an initial power-law increase according to [Of 

I« « 0.78 7 - 4 n 1/3 , 7 < 1 , 1< n < 7 3 . (5.7) 

At n « 0.297 3 it reaches a maximum and then decays exponentially according to 



IV - ~~~~~t exp f-|^) , 7 » 1 , n » 7 3 • (5-8) 

This behavior is illustrated in Fig. |5](c) for the case s = 10 (full line). The analytical results 
Q5.7| ) and ( |5.8| ) (dashed lines) are also shown in Fig. |5](c). They compares well with the data 
in the appropriate limits. We now show that the behavior of in conjunction with the 
behavior of g n leads to substantial suppression of synchrotron radiation for cold beams. 
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FIG. 5. The normalized partial powers In as a function of n for (a) s = 0.1, (b) s = 1 (b) 
and (c) s = 10. The analytical results (|5Jj), (|5.6|) , ( |5.7D and (|5.8|) are also shown in the respective 
panels. 

We first discuss the case of a coasting crystallized linear chain at T = 0. It consists of 
N equi-spaced particles according to 9j = 2nj/N, j = 1, 2, . . . , N. For g n we have the result 
( p.20| ). For the suppression factor a we obtain in this case 
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a(N,(3)=N £ l£ N . (5.9) 

m=l 

We saw above that independently of s the normalized partial powers 1^ always decay 
exponentially for large enough n. Thus, there is always an Nq such that 

a^NI N (5.10) 

to a very good approximation. Therefore a is exponentially small for N > Nq. In other 
words: For large enough particle number we obtain exponential suppression of synchrotron 
radiation independently of the relativistic regime of the beam. This result is illustrated in 
Fig. ^ (r = oo case). It shows the suppression factor for s = 0.1, 1 and 10 as a function 
of the particle number N. In all three cases we indeed obtain exponential suppression as 
predicted from the structure of (]5.9|). 
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FIG. 6. Suppression factors for the crystallized chain for three different plasma parameters 
(r = 10,100,oo) in the three relativistic regimes: (a) s = 0.1, (b) s = 1 and (c) s = 10. The 
theoretical curves correspond to equations ( 5.11[ ), ( |5.12| ) and ( |5.13| ). 



Using Eq. fl5.10|) and the above expressions for I n in the relevant relativistic regimes, we 
obtain explicit analytical formulae for a(N,/3): 



a(N,{3) 



3iV§ 



2iV 



4-1 , s < 1, A^> 1 , 



(5.11) 
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3 / i \ 2JV 

a(N,f})a ^ , ^ , s<l,iV» 7 3 , (5.12) 

a(iV,/3)» 9 expf-^ , s » 1 , N » 7 3 • (5.13) 

Fig. ^| shows that the analytical formulae are very good approximations of the numerical 
data in their respective ranges of validity. 

Next we consider the linear chain at finite temperature. In this case the form factor 
( j3.19| ) applies. Because of the structure of (|3.19|) and the asymptotic exponential decay of 
1^ for large n we can compute the asymptotic behavior of a(N, (3) for large N. Using ( |3.22p 
we obtain 

-j OO -I OO 1 

« W ft = m E ^ w N 9 ^ E ^ = ^7^ « — , AT » 7 3 . (5.14) 

iV n=l n=l /X 

Thus, for large iV and in all three relativistic regimes, the asymptotic suppression is inde- 
pendent of N and saturates at a = 1/(2T). This behavior is illustrated in Fig. |6] which shows 
the suppression factor for T = 10, 100 and oo for all three values of s considered. The onset 
of saturation in the vicinity of some N = N c is physically clear because of the following 
reason. Finite T corresponds to a finite temperature which furthermore corresponds to a 
finite correlation length of the particles in the linear chain. But since the suppression of the 
synchrotron radiation is a coherent process it is intuitively clear that no further suppression 
can be achieved once the total particle number exceeds the correlation length. Consequently 
the suppression effect has to saturate. 



In Sec. |nj we pointed out that measuring the depth of the correlation hole in g" y for 
small values of n defines an experimental method for measuring the plasma parameter of 
the beam. Since the saturation value of a depends only on T, measuring the suppression 
factor for large N defines yet another experimental procedure for measuring V . 

The existence of a finite correlation length at finite temperature provides an argument 
why it is not necessary to maintain coherence over the whole circumference of the storage ring 
in order to observe the suppression effect. It is enough to work with bunches whose length 
is smaller or of the order of iV c to observe suppression of synchrotron radiation. From the 
mathematical point of view this is also evident since we saw that for finite temperature only 
lower harmonics of g n are affected by the bunching, hence gm y N nnch ~ 9rwNi w = 1, 2, ... for 
large enough N, and therefore we expect similar suppression as for the coasting case. This 
result is very important for practical applications of the suppression effect. It means that the 
coherence does not have to be maintained over the whole extent of the ring, which sometimes 
can amount to hundreds of meters and more. It is enough to maintain the crystalline 
structure over small angular distances (bunches) in order to exploit the suppression effect 
in possible technical applications. 

It is also clear by inspection of Fig. || and of Fig. |2](b) that for large enough N substantial 
suppression of synchrotron radiation can be achieved for liquid beams. This, again, is 
important since modern electron coolers are close to providing a liquid beam of electrons. 
Thus it may soon be possible to check our theory with the help of liquid electron beams. 
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VI. DISCUSSION, SUMMARY AND CONCLUSIONS 



The suppression of the radiation of geometrically ordered charges was first noticed by 



J. J. Thomson ||20|| . He employed this effect for motivating the stability of atoms, which, 
according to classical theory, should radiate and decay. Suppression of synchrotron radiation 
in the context of accelerators was first noted by L. I. Schiff ||12|| . But in Schiff's time a 



mechanism for establishing the order in a beam of charged particles was not available. Only 
recently, with progress in the cooling of beams by electrons and lasers is it possible to 
envision the production of crystallized beams whose synchrotron radiation is exponentially 
small. It should be born in mind, however, that synchrotron radiation is not very important 
for heavy ion beams that can easily be cooled with electrons and lasers. Dramatic effects are 
expected to occur only for crystallized electrons where the synchrotron radiation is orders of 
magnitude stronger. The draw-back is that electrons cannot be cooled directly with lasers. 
We hope, however, that this paper will stimulate experimentalists to develop cooling schemes 
for electron beams. One possibility would be to use sympathetic cooling of electrons with a 
beam of heavy ions that can be cooled by lasers. 

The paper discusses various forms of ordered beams that may occur in practice: Gaseous, 
liquid and crystalline, coasting and bunched. It is pointed out that the suppression effect 
occurs on two levels: In the form factor of the beam and in the total radiated power. While 
the modifications in the form factor may be used as a diagnostic tool for inferring the 
thermodynamic state of the beam, the suppression of the total power may eventually lead 
to the construction of small-sized cyclic electron accelerators. 
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